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. The X-ray cascade from antiprotonic atoms was studied for 208 Pb and 209 Bi. Widths and shifts 

' °f the levels due to the strong interaction were determined. Using modern antiproton-nucleus 

optical potentials the neutron densities in the nuclear periphery were deduced. Assuming two 
£SJ . parameter Fermi distributions (2pF) describing the proton and neutron densities the neutron rms 

radii were deduced for both nuclei. The difference of neutron and proton rms radii Ar np equal 
to 0.16 ± (0.02) stat ± (0.04) S j, s < fm for 208 Pb and 0.14 ± (0.04) stat ± (0.04) SJ/St fm for 209 Bi were 
determined and the assigned systematic errors are discussed. The Ar np values and the deduced 
shapes of the neutron distributions are compared with mean field model calculations. 
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I. INTRODUCTION 



Beginning more than ten years ago, we have performed an experimental study of the medium-heavy and heavy 
antiprotonic atoms using the slow antiproton beam from the Low Energy Antiproton Ring (LEAR) at CERN. The 
main objective of our program was to obtain information on the neutron distribution at the nuclear periphery and to 
provide data useful for deducing the antiproton-nucleus optical potential parameters. 

Two experimental methods were employed. First, using the so called "radiochemical method" we have investi- 
gated [l|, 0j S 0] the ratios of peripheral neutron to proton densities at distances around 2.5 fm larger than the nuclear 
charge half-density radius || . The method consisted in measuring the yield of radioactive nuclei having one proton 
or one neutron less than the target nucleus, produced after antiproton capture, cascade and annihilation in the target 
antiprotonic atom. The experiment yielded 19 density ratios (proportional to the "halo factor" /h a io defined below), 
which were subsequently employed to deduce the shape of the peripheral neutron distribution. 

The second method consisted in measurements of the antiprotonic-atom level widths and shifts due to the strong 
antiproton-nucleus interaction. These observables are sensitive to the interaction potential which contains, in its 
simplest form, a term depending on the sum of the neutron and proton densities. The level widths and in a number 
of cases also the level shifts were measured for 34 antiprotonic atoms (in some cases for different isotopes of the same 
element). 
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The rich harvest of the two methods employed which are sensitive to the neutron and proton density ratio and the 
sum of these densities, has allowed to derive a number of systematic conclusions on the nuclear periphery properties 
presented in a series of summary and analysis publications [1, 0, H, Moreover, our data were used to determine 
the antiproton-nucleus optical model parameters through global fits of p X-rays and halo factors [l(J HH with a 
substantially larger and more precise database than employed in previous approaches [l2l . [l3j . 

Besides our summary papers, after the end of the antiprotonic X-ray (PS209) experiment, we prepared more detailed 
reports, containing information on experimental procedures and their analysis for some cases studied [Til. [TEl HtH [l7T] . 
The present article, dealing with 208 Pb and 209 Bi antiprotonic atoms, is the next in this series. 

During the last years it was shown that properties of the 208 Pb neutron distribution can be correlated with a 
number of quantities in various fields. In particular the knowledge of the difference Ar np between the rms radii of 
neutrons and protons in this nucleus constrains the symmetry energy of nuclear matter and therefore is reflected in 
the neutron Equation of State (EOS) [3, 0, HO, HH, H3] • The neutron EOS models, in turn, are used to calculate 
the properties of neutron stars, such as their radii and proton fraction [HI, [Hj]. However, not only the first moments 
of the neutron density distributions, but also their shapes are of considerable interest, e.g. in the determination of the 
isovector potential parameter of the pion-nucleus s-wave interaction in nuclear matter [261 ] or in the calculation of the 
lepton flavor violating muon-electron conversion rate [271 ]. There is also a certain dependence on the radial neutron 
distribution in the proposed determination of the 208 Pb neutron rms radius through parity violating electron-nucleus 
scattering 

Experimentally, the Ar„ p value in 208 Pb was previously determined u sing hadronic probes (elastic scattering, 
inelastic scattering exciting GDR), reported in Refs. [29], |3Cj, HH, [H, [H, [34], [3f| [36| and discussed in There were 
also some attempts to deduce higher moments of the neutron distribution from the hadron scattering experiments [13] 
(see also [36|). On the other hand, the abilit y o f the medium-energy elastic proton scattering data to determine the 
neutron distributions was recently contested [381 ]. 

The measurement in 209 Bi offers other advantages. It is the only experiment that allows to see an even-odd isotopic 
effect in heavier nuclei, in this case due to the loosely bound proton in 209 Bi. One difficulty in the way of analysis is 
the calculation of the hyperfine-structure that is comparable to the strong interaction shift and broadening. After this 
is done, it turns out that the level shift in 208 Pb is repulsive (as most of the lower shifts), but the level shift in 209 Bi 
is attractive. This finding is open to interpretation. Here we pursue the view that it is related to a pN quasi-bound 
state, which is important in cases of loosely bound valence nucleons. 

II. EXPERIMENTAL METHODS 

The heavy antiprotonic atoms 208 Pb and 209 Bi were investigated during the experiment PS209 at CERN in 1996 
using antiprotons of momentum 106 MeV/c. Table U gives the target properties and the number of antiprotons used 
for each target. 

The antiprotonic X-rays emitted during the antiproton cascade were measured by three high-purity germanium 
(HPGc) detectors. Two detectors were coaxial with an active diameter of 49 mm and a length of 50 mm (relative 
efficiency about 19% and 17%, respectively), and the third one was planar with 36 mm diameter and a thickness of 
14 mm. The detectors were placed at distances of about 50 cm from the target at angles of 13°, 35° and 49° towards 
the beam axis, respectively. The detector-target distance was adjusted to obtain a good signal-to-noise ratio and to 
decrease at the same time the background produced by pions from the annihilation processes. 

More details concerning the experimental methods, the detector calibration and the data reduction may be found 
in our previous publications [14j . Il5l . Il6l . Il7| . 

III. EXPERIMENTAL RESULTS 

The strong interaction between antiproton and nucleus causes a sizeable change of the energy of the last X-ray 
transition from its purely electromagnetic value. The nuclear absorption reduces the lifetime of the lowest accessible 
atomic state (the "lower level" , which for lead is the (n, I = 9, 8) state) and hence this X-ray line is broadened. 
Nuclear absorption also occurs from the next higher level ("upper level") although the effect on level energy and 
width is generally too small to be directly measured. The width of the (n, I = 10, 9) level was deduced indirectly by 
measuring the intensity loss of the final X-ray transitions. The level scheme for the antiprotonic Pb atom with the 
observables of the X-ray experiment is shown in Fig. [TJ 

The X-ray spectrum measured with antiprotons stopped in 208 Pb is shown in Fig. [2j Those lines in the spectra 
that are not broadened were fitted with Gaussian profiles. The lowest observable LS-split doublet lines (n = 10 — > 9), 
which are significantly broadened, were fitted with two Lorentzians convoluted with Gaussians (Fig. [3]). 
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The measured relative intensities of the antiprotonic X-rays observed in the investigated lead and bismuth targets 
are given in Table QT] These intensities were used to determine the feeding of the consecutive n levels along the 
antiprotonic-atom cascade. This is shown for 208 Pb in Fig. [4] 

Table Hm gives the measured shifts e, defined by e = E em — E exp , where E cxp the experimental value for the transition 
energy and E em the energy calculated without strong interaction [3i|. For the lead (n, I = 9,8) levels the shifts are 
clearly repulsive, whereas for bismuth the levels shifts are consistent, within the errors, with zero (Fig. [5|). However, 
as discussed in the next sections, there exists an additional repulsive shift due to the hyper fine structure. For the 
(n, I = 10, 9) levels in lead the shifts are repulsive, as in the case of the (n, I = 9, 8) levels, but the shifts are smaller. 
Tables [IV] and |V] give the measured widths. As indicated above, the widths of the (n, I = 10, 9) levels were derived 
from the intensity balance of transitions feeding and depopulating these levels. Contributions of parallel transitions to 
the measured intensities were obtained from cascade calculations (see [HI). The rates for radiative dipole transitions 
were calculated with the formulae given in Ref. (40| . The Auger rates were derived from the radiative rates and from 
cross sections for photoeffect using Ferrell's formula [4l|. The width of the levels (n, I — 10,9) are larger for 209 Bi 
than for 208 Pb (Fig. [5j . This is due to the hyperfine structure in Bi, which will be discussed below. 

IV. ANALYSIS AND DISCUSSION OF THE 208 Pb MEASUREMENT 
A. Introductory statements 

The analysis of the presented antiprotonic 208 Pb atom data is based on some assumptions that are briefly mentioned 
here and discussed below. We assume and show by a comparison with model calculations that charge, proton 
and neutron distributions in this nucleus can be well approximated by two-parameter Fermi (2pF) distributions: 
p(r) = po • {l+exp(^ :: ^)} _1 , where c is the half-density radius, a is the diffuseness parameter, and po is a normalization 
factor. In particular, calculating the neutron rms radius from the antiprotonic X-ray data sensitive to densities at 
distances around 1.5 fm larger than the half-density charge radius @, we extrapolate the experimental density well 
into the interior of the nucleus. As shown below (Sec. II V C[) . this assumption is in reasonable agreement with the 
density shapes calculated in terms of the mean field models. 

In evaluating the observables of the antiproton-nucleus interaction the important question of the ratio of the 
annihilation probability on a neutron to the one on a proton arises. In the simplest nuclear optical potentials this ratio 
is given by the ratio of the imaginary parts of the effective scattering lengths, R = Imb^ / Imb^ . The experimental 
determination of this quantity [42|, |43| gave R = 0.63. 

In spite of this observation, a value R=l was assumed in the optical potentials proposed in Refs. 
Our analysis of the antiprotonic X-ray data, together with the radiochemical experiment, also indicated 
consistency of these two methods when R — 1 was chosen. Therefore this value is also adopted in the present data 
evaluation, as discussed in the following sections. 

B. Charge and proton distributions 

It is generally assumed (see e.g. Ref. (45|) that the charge rms values are known for the stable nuclei with remarkable 
precision, about 0.3%. The same belief is often projected on the charge distributions. In Fig. [6] we show (as it was 
already observed in [l7j ) that this is not the case for the radial distances where the antiproton-nucleus interaction 
takes place in 208 Pb (about 7 to 10 fm away from the nuclear center, see below). In this figure the charge density of 
208 Pb, tabulated in a number of compilations, is compared with the most recent one by Fricke et al. (45|. Neglecting 
the oldest tabulation, differences of up to 50% are observed between Fricke et al, Jager et al. [46| and de Vries et 
al. [47J for radial distances close to 10 fm. We consequently use the Fricke charge distribution in this work. 

Experiments using electromagnetically interacting probes give charge density distributions or rms charge radius 
values (e.g. l45l I471B whereas point proton distributions are needed when Batty 's zero-range antiproton-nucleus op- 
tical potential |12ll is used for the analysis of the experimental data. For the finite-range version of the p-nucleus 
potential jHJ, [Il[ these point distributions are folded over an interaction range. 

In Ref. [48| the analytical formulae to transform the 2pF charge distribution to the 2pF point distributions of the 
proton centers were presented. We have previously used them in our data analysis, presented e.g. in Ref. Q. Similar 
analytical formulae were recently given in Ref. [4!|. In order to transform the 208 Pb proton charge distribution of 

Ref. [ill, we have used the proton charge rms radius < r 2 > = 0.875 fm obtaining 2pF point proton parameters 

of c p — 6.684 fm, a p = 0.446 fm and rms r p = 5.436 fm. 
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C. Calculated mean field neutron and proton distributions 



The proton and neutron distributions in the doubly magic 208 Pb nucleus were subject of a large number of theoretical 
investigations. In this paper we select Skyrme-Hartree-Fock (HF) and the Hartree-Fock-Bogoliubov models, namely 
those with the SkP (HFB) and SkX (HF) [52| parametrization, both reproducing the 208 Pb charge (proton) radius 
and neutron binding energy remarkably well. (It has, however, recently been shown that the SkP Skyrme model may 
diverge for some nuclei if calculated to sufficient accuracy [53J| ) . The third self-consistent mean filed model considered 
here belongs to the framework of the relativistic mean field theory (RMF) with the recent DD-ME2 parametrization 
of the effective interaction (54|. Although in fitting the DD-ME2 parameters the 208 Pb Ar np value (of 0.20 fm) was 
used to adjust the interaction parameters, the shape of the neutron distribution was obtained from the calculation. 

Figure [7j left panel, presents the proton and neutron distributions of 208 Pb as calculated using the DD-ME2 
parametrization. This and two other distributions were approximated by 2pF distributions fitted to the theoretical 
densities. Satisfactory fits were achieved for radii between 1 and 10 fm (local differences between theoretical and 
fitted densities were less than 10% for protons and less than 4% for neutrons in the radius range 2-9 fm, Fig. right 
panel). Figure [5] shows the summed neutron and proton densities for the three forces considered. The relationship 
between the neutron densities, rms radius and equation of state is being investigated with a larger set of mean-field 
models in [55| . 

Table I VII gives the results of the fitting procedure. The rms radii calculated using a and c values of the 2pF are 
close to those obtained using the theoretical distributions directly. The neutron distributions are close to the "halo 
type" 0, with Ac np = 0.02 fm for SkP, Ac np = 0.07 fm for SkX and Ac„ p =0.05 fm for DD-ME2, respectively. 
This is illustrated in the left-hand part of Fig. [51 showing the normalized neutron to proton density ratio obtained 
from the density distributions of the discussed models. The figure presents also the "pure halo" distribution with 
Ar nD =0.16 fm and Ac nlD =0 fm. 



D. Antiproton-nucleus optical potentials 

The standard potential for hadronic atoms [l3| are composed of two terms 

V opt = V s (r) + Wp(r)V (1) 

involving proton and neutron components. Both terms, the local Vs and the gradient Vp are expected to have the 
folded form 

V s , P (r) = -—b s , P / du/(u)p(r - u), (2) 

which involves nucleon densities p, folded with a (usually Gaussian) function / of some rms radius ,r range , of the 
range, "effective lengths " 6, and some weak dependence of the reduced mass p N jj on nuclear recoil. It turned out 
already in a first analysis for heavy atoms that an independent determination of b s and b P is not realistic. A 
simplified result in Ref. [l3[ gave bs — —2.5(3) — i 3.4(0.3) fm with zero r range , which we call the Batty potential 
and use in our calculations. The recent phenomenological best fit (Friedman potential) is obtained with a single 
term of bs (= bo) = —1.3(1) — £1.9(1) fm for both neutrons and protons and r range — 1.04 fm [l(j. Within all 
these calculations the X-ray data suggested no significant differences in the the values of the pp and np lengths. No 
relation of the phenomenological 65 to the NN scattering parameters has been established. Another recent analysis 
[liT | attempts to find such a relation from the analysis of the lightest atoms H, D, He and of scattering data described 
in terms of the Paris NN potential. One important difference arises in equation (|2|), which also contains nucleon 
recoil terms not required in a phenomenological approach. The nucleon recoil term constitutes about a quarter of the 
total potential and depends on the state of the nucleon. This leads to a more complicated parametrization, which 
will not be repeated here. In this potential [llj one obtains roughly 65 = —1.7 — i 0.9 fm, bp = — z 0.4 fm 3 with 
Trange = 0.8 fm and the absorptive part of these parameters compares well with the pp and ph scattering data. 

Presently available optical potentials are unable to reproduce the level shifts in Pb. This reflects a more general 
difficulty related to uncertainties of the real part of V opt . More specific difficulties such as pN-quasi-bound states or 
long-range pion exchange forces in the pN system have already been discussed in Refs. [Til. f56| . 



E. Annihilation probability 



The probability P n ,i{r) of the nuclear capture from a given atomic state (n, I) is defined (for local optical potentials) 
as P n ,i{f) =| 4>n.i(r) | 2 (ImV opt ) ■ r 2 , where 4>n,i is the antiprotonic wave function and r is the radial distance from 
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the nuclear center. 

The most probable value of this probability was calculated for some cases in Ref. [5j and the results of these 
calculations were confirmed by the analysis given in Ref. [1(3]. The calculations of the P n ,i(r) distributions for 
208 Pb were preformed using Batty potential and density distributions discussed in Sec. IIV CI Figure [TU] shows this 
annihilation probability distributions for DD-ME2 densities and Tabic Ivm gives parameters of this distribution. Wider 
distributions (by about 60% for upper level and 20% for lower level) with more pronounced tails for larger radii are 
obtained using finite range optical potentials. 

F. Interpolated halo factor 

In our previous analysis of the antiprotonic X-ray and radiochemical data the crucial information used in deducing 
the shape of the neutron distribution ( "halo type" or "skin type" , see 0] for the definition) was the experimentally 

determined halo factor /h a io = y^l-i) ' lt t ' ^' nere ^ are t ^ ie J^ds f° r tne — 1 nuclei, Z t , N t and A t are the 
target proton, neutron and mass numbers, respectively, and R was defined and discussed above. We have shown 
previously that in a number of cases for which this halo factor could be measured, it indicated that the corresponding 
2pF neutron distribution is close to the "halo type" , i.e. with equal proton and neutron radius parameter (c n =c p ) and 
larger diffuscness parameter for neutrons (a„>a p ). Again, this observation is in fair agreement with the mean-field 
calculations of the nuclear densities (see Sec. IIV C|) . 

In the case of the 208 Pb nucleus, the experimental determination of the halo factor by the radiochemical method 
was not possible as one of the (At — 1) isotopes ( 207 Pb) is not radioactive. In order to have at least some indication of 
its value the 208 Pb halo factor was deduced by interpolation between /halo values for other nuclei, which are plotted 
either as a function of the neutron binding energy [2| (see Fig. or of the asymmetry parameter 5 = (N — Z)/Z 
[571 ]. (Note that, as discussed in Ref. 0, the halo factors published by us before this reference should be multiplied 
by 0.63). The interpolated / halo for 208 Pb found in this way is 2.8 ± 0.4. 

This interpolated value can be compared with the results of Bugg et al. [42|, where the idea of the neutron halo 
factor was introduced for the first time to interpret the ratio of charged pions generated by antiproton annihilation 
in various targets, including 208 Pb. In that work the halo factor, defined as fH alo = jv(p^) W ^ ( wnere N(np) and 
N(j>p) are the number of p annihilations on peripheral neutrons and protons, respectiv ely) was measured for the 
208 Pb nucleus by detecting the charged pions emitted after antiproton capture in nuclei [421 ]. The result obtained, 
fhaio = 2.34(50), is based on the assumption that R = 0.63, with this latter number extracted from the p capture 
in carbon. The result has been subject to some criticism for neglecting the final state interactions and a possible 
dependence of R on the actual nucleus. The more recent Obelix experiments determined R = 0.48(3) from low- 
energy capture in He [58j, while the best fit optical potential requires R « 1 [10]. This discrepancy is resolved if one 
realizes that capture in He involves mostly S-waves and capture from high angular momentum states in heavier nuclei 
involves mostly P (or higher) waves in the NN system [59| . We refer to Ref. [5j| for some details of the calculation, 
which estimates R « 0.9 — 1.0 in the lead region and for a new analysis of the Bugg result. R = 1 is used in these 
calculations. This yields /£ a ; (Pb) = 1.8(4). One can obtain the average radius of the absorption r 77 via calculations 
of the final state pion interactions and comparison of the final experimental and calculated pion spectra. This yields 
r T = c c h + 1.35 fm. 

G. Experimental results analyzed by optical potentials 

1. The zero NN range antiproton-nucleus potential of Batty 

Our previous results [3] for the neutron distribution in 208 Pb were obtained using a zero-range NN force Batty 
potential ^2Tll3|. This potential, obtained by a fit to the antiprotonic X-ray data mainly on light nuclei was determined 
before our results were published. It has, therefore, the evident philosophical advantage when compared to the recently 
published potentials strongly [11] or completely [l(| relying on our results, including the 208 Pb nucleus. 

The neutron-distribution parameters deduced using Batty's potential are shown in Table IVIl Calculations were 
done with a code based on the work by Leon [60] . The Ar np is calculated as the difference between the rms radii of 
the corresponding 2pF point proton (Fricke) and point neutron distributions under the assumption of a "halo type" 
distribution (Ac np =0). The Ar„ p value of 0.16 fm differs by 0.01 fm from the previously published one Q; this is due 
to the updated value [Hot ] of the electromagnetic proton rms radius used in the transformation from charge to proton 
densities. 
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As indicated above, our experimental data yield Ar np — 0.16 fm under the assumption of c n = c p . In Figs. ll2l and ll3l 
we show how the relaxation of this condition would influence the difference of the Fermi-distribution parameters and 
the rms radii difference for the 208 Pb nucleus. 

The data in Fig. fTS] present the change of Aa np when we allow the Ac np to change while still being in agreement 
with the experimental level widths. It is seen that with the extreme neutron-skin assumption (identical proton and 
neutron diffuseness, Aa np =0) the neutron-proton difference of the half density radii Ac np should be close to 0.8 fm. 
As shown in Fig. [T3] such a large value of this difference would lead to Ar np close to 0.6 fm. 

In Ref. Q we have discussed the results for Ar np in 208 Pb obtained by using hadron scattering data. The weighted 
average of six experimental results, obtained between 1979 and 2003, is Ar np = 0.16 ± 0.02 fm @. It is in excellent 
agreement with the result obtained from the antiprotonic X-rays with the zero interaction range Batty potential. The 
gray band in Fig. [13] indicates the error margin of the weighted average of the hadron scattering experiments, allowing 
a difference in the Ac np value between the 2pF distribution of neutrons and protons in 208 Pb to be 0.08 fm at most. 

The upper limit for Ac np can also be estimated comparing the calculated neutron to proton density ratio with the 
interpolated halo factor of 208 Pb. This density ratio is plotted in the right panel of Fig. [9] as a function of the radius 
for Ar„ p =0.16 fm (Batty potential X-ray value for Ac np =0 and average value from the hadron scattering data). Two 
halo factors are shown in this figure: one resulting from the pion emission experiment and the other one from the 
interpolation (see Sec. IIVF|) . Although the pion experiment is not limiting the Ac np value, the interpolated /halo 
clearly indicates that Ac np has to be smaller than 0.1 fm. This determines the systematic error of the Ar np value to 
be equal to 0.04 fm. If Ref. [46[ instead of Ref. [45| is taken for the charge distribution, Ar np =0.12 fm is obtained, 
i.e. the systematic error is also 0.04 fm for this case (cf. also results of Ref. [llj in Sec. IIV G~3)) . 

The assigned statistical and systematic error for the Ar np value indicates about 1% uncertainty in the determination 
of the neutron rms radius in 208 Pb from the antiprotonic atom data. This value is comparable to the expected precision 
of the parity violation measurements 28] of this quantity. 

The comparison of the experimentally determined level widths and shift with the theoretical proton and neutron 
distributions (see Sec. IIV Cjl using Batty's potential is shown in Table IVIIII The level widths calculated with the 
SkP and SkX distributions are too small (by 12% and 27%, respectively) whereas the DD-ME2 distributions result 
in widths close to the experimental values. As indicated in Sec IIVD1 the value of the shift is not reproduced for any 
theoretical distribution. 

It is interesting to note that the SkX interaction with the value closest to the experimental Ar np leads to level 
widths which are clearly below the observed ones. It indicates that for this interaction the nucleon density decreases 
too fast with r. This was illustrated in Fig. [5] by a comparison of the summed neutron and proton densities for the 
three forces considered. As previously discussed in [36| . too small proton but especially neutron diffuseness exhibited 
by the SkX model is the reason for this behavior (cf. Tabic Wl} . Our antiprotonic atom data presented in Table IvnTl 
are therefore a confirmation of the conclusions drawn from the analysis of nucleon elastic scattering, presented in 
Ref. [£|. 

2. The finite NN range potential of Friedman 

Contrary to the Batty potential [HI, EH , the finite-range antiproton- nucleus interaction potential recently proposed 
by Friedman et al. [lfj was based almost completely on the PS209 experimental data, including the antiprotonic atom 
level widths and shifts of 208 Pb and 209 Bi reported in this publication. 

In order to obtain a similar relationship as that shown for Batty's potential in Fig. [121 the 2pF distribution of 
protons (see Table fVT]) was folded over the interaction range with 1.04 fm rms radius. The (folded) neutron density 
was fitted using the optical potential parameters from [l(| to the experimental level widths of 208 Pb, varying Ac^ p . 
The Aa^ values of the folded densities were deduced from the fit and are shown in Fig. [T3J Figure [151 similar to 
Fig. [T31 shows the resulting Ar^ p values of the folded density distributions as a function of the folded Ac^ p values. 

It can be shown that the transformation from the point-like nucleon distributions to the folded ones increases the 
Aa^p values by about 15%, leaving Ac np and Ar np approximately unchanged. Therefore, without any supplementary 
conditions these results indicate that using the finite-range version of the optical potential as proposed by |l0j . the 
Ar/ r »Ar r value obtained for 208 Pb is between 0.11 fm (Ac f np =0 fm) and 0.38 fm (Ac£ p =0.54 fm). 

In Ref. 10] it was shown that if all antiprotonic data are presented in the form Ar np = a(N — Z) / A + [3 Q the 
global fit to these data allows a values between 0.9 fm and 1.3 fm (our previous data in 0] analyzed in terms of 
the point proton optical potential gave a = 0.90 ± 0.15 fm). This a range gives Ar np values between 0.154 fm and 
0.240 fm. 

Another limit of the Ar np values could be obtained with the help of the weighted average of the hadron scattering 
data, giving Ar np =0.16±0.02 fm Q. Two standard deviations of this average (cf. Fig. fTS"]) put the lower and upper 
limit of the Ar nv as 0.12 fm and 0.20 fm. 
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Taking the lower allowed value from the global fit and the upper one from the hadron scattering gives 
Ar£ p «Ar nj) =0.17 fm as the result of the finite range Friedman potential with the estimated error of ±0.02 fm. 

This result implies Ac£ p ~Ac rap of 0.13 fm, i.e. the neutron distribution essentially of the "halo type" but with a 
small contribution of the "skin type" . Such a distribution is shown in Fig. [9] by the curve labeled "G" , which lies 
slightly below the lower limit of the interpolated halo factor. 

The comparison of the experimentally determined level widths and shift with the theoretical proton and neutron 
distributions using the Friedman potential is given in the lowest part of Table IVIIII The calculated level widths have 
about 5% uncertainty due to the folding procedure applied. Within these errors calculated level widths are identical 
to those obtained using Batty's potential. 



3. A constraint finite NN range potential 

A parallel study with the potential from Ref. 11] generates two solutions Ar„ p =0.16(3) fin for the electron scattering 
charge [46[ and/or Ar„ p =0.22(3) fm for the muonic charge density (45|. These solutions favor halo type neutron 
densities, both are characterized by large \ 2 values due to the poorly reproduced level shift. 



V. ANALYSIS AND DISCUSSION OF 209 Bi DATA 



Since the 209 Bi nucleus has the spin 1 = 9/2 and magnetic moment /i = 4.08 • /in (where /ijy is the nuclear 
magneton), the antiprotonic atom levels are split. The related hyperfine shifts are smaller than the fine structure 
(f.s.) splitting but comparable to the strong interaction shifts. The standard formula [6l{ . extended to the case of an 
anomalous magnetic moment (62j . gives 

F a 9p9Bi [F(F + 1) - J(J + 1) - 1(1 + 1)] 3 

EhSs = ~limM~ J(J+l) n 3( Z + l/2) imZa) ' (3) 

where M is the proton mass, m is the reduced mass, gl = \i in nuclear magnetons. The last transition in Bi is split 
into 10 dominant components in the upper f.s. state and into 9 components in the lower f.s. components. These 
correspond to different values of the total spin of the system F = J + I, .... J — I. Assuming a statistical population 
~ (2F +1), the observed spectral line becomes asymmetrical. One obtains an overall 50 eV repulsive shift of the 
centroid and an additional "broadening" of 150 eV. This yields a lower width of 350(50) eV and an attractive lower 
shift of 37(53) eV, generated by strong interaction. (The upper level width averaged over the fine structure components 
(see Table fV| is equal to 6.9(1.3) eV). One thus faces a sizable isotopic effect between attraction in 209 Bi and repulsion 
in 208 Pb atoms. The difference is related to the weakly bound valence proton in this Bi isotope. As discussed on 
previous occasions, there are indications of a quasi-bound state in the pp system just below the threshold. Such a 
state generates few distinct cases of an anomalous behavior of level shifts in nuclei with loosely bound nucleons . 
For a more quantitative discussion of these phenomena we refer to a parallel publication [ll[ . 

Assuming (as for 208 Pb) Ac„ p =0 fm, the difference of the neutron and proton diffuscness parameter was fitted to 
level widths using Batty's potential. The Ar np value obtained is equal to 0.14±0.04 fm, lower by 0.04 fm than the 
previously published one Q . In Ref. [7| the hyperfine splitting of the lower level discussed above was not taken into 
account. 



VI. SUMMARY AND CONCLUSIONS 



In this article we have presented an analysis of the nuclear structure information extracted from the studies of 
the antiprotonic atoms p- 208 Pb and p- 209 Bi. The experimentally determined level widths and shifts of these atoms 
at the end of the antiprotonic cascade depend on the antiproton-nucleus interaction potential. In turn, the crucial 
ingredient of this potential is the nucleon density at the radial distance where the antiproton-nucleus interaction 
occurs. Therefore, as it was shown already in our first publications in this series P, the study of antiprotonic 
atoms may constitute a powerful tool for the extraction of information on the properties of the nuclear periphery. 

In the analysis of the antiprotonic atom data presented here we have been essentially using two optical antiproton- 
nucleus potentials, proposed by Batty et al. (l2l . [l3| and Friedman et al. flo| . respectively. The Batty potential, 
now more than ten years old, was obtained by fitting the potential parameters to the 33 level widths and 15 level 
shifts of antiprotonic atoms published at that time, mainly of light and a few intermediate-mass nuclei. The fits 
were performed with a zero-range antiproton-nucleon interaction. Although this unphysical assumption is presently 
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avoided jlOL lllj as leading to worse fits than the finite range potentials, we still pursue our data analysis using the 
Batty prescription in its simplest form. Our arguments are that although "point nucleon distributions" and "zero 
range interaction potentials" are probably an oversimplification, the obtained parameters are deduced from the fit to 
the experimental data. It may be expected that this fact somehow in an automatic way introduces the corrections 
of the method deficiencies. Moreover, as already mentioned at the beginning of Sec. IIV G 1[ the Batty-potential 
parameters were obtained before our antiprotonic atom data were available, ensuring the interpretation of the results 
to be independent of the interpretation tools. 

Recently, an antiproton- nucleus optical potential with finite interaction range was proposed by Friedman et al. [Ic| . 
It was shown that the 90 data points from our PS209 X-ray experiments, together with 17 data points from the 
radiochemical experiment, determine an attractive and absorptive p-nuclcar isoscalar potential, which fits the data 
well. 

However, as the 208 Pb antiprotonic X-ray data were used in the determination of the Friedman potential, we have 
tried to show in the analysis of the experiment with this potential what can be deduced on the neutron distribution of 
this nucleus on a more general ground. To this end we have used the information on the trend of the Ar„ p values as 
a function of the asymmetry parameters (N — Z)/A, allowed by the potential of Ref. [lfj, and previously analyzed Q 
results of the hadron scattering experiments. 

Another fit with a finite-range potential was recently proposed (Ref. [ll|, see also Sec. IIV G 3|) . The reader is 
referred to the original reference for the discussion of the global fit to the antiprotonic X-ray data performed and the 
Ar np values deduced for 208 Pb. These values are used in the present publication to estimate the systematic errors. 

It was shown in this paper that by an analysis based strictly on the experimental antiprotonic level widths one 
would be unable to propose meaningful limits for the Ar np value in 208 Pb. Applying the Batty potential with a pure 
"halo shape" of the neutron distribution (Ac np =0 fm) leads to Ar np =0.16 fm. The lowest limit of the interpolated 
halo factor would allow at most Ac np =0.1 fm, i.e. Ar np =0.20 fm, similarly to the 2a uncertainty of the average Ar np 
value deduced from the hadron scattering experiments. The analyzed theoretical proton and neutron distributions 
using HF, HFB and RMF models give a Ac np value of 0.07 fm at most. If the charge distribution from Ref. [46| is 
used instead of that from Ref. [45|, Ar np =0.12 fm is obtained. We conclude that our experiment interpreted using 
Batty's potential with the supplementary information given above leads to Ar„ p =0.16 ± (0.02) stat ± (0.04) sast fm in 
the 208 Pb nucleus. A value for Ar rap that is only 0.01 fm larger is deduced from the analysis using the Friedman 
potential. 

Significant results were obtained applying the Batty and Friedman potentials with the theoretical proton and 
neutron distributions to get the antiprotonic 208 Pb atom level widths and shift. As it was discussed in Sec. IIV Dl the 
experimental level shift was not reproduced. On the other hand for each theoretical distribution both potentials give 
almost identical widths in spite of the fact that one is "zero range" and the other one is "finite range" interaction 
potential. The calculated widths are smaller in some cases than the experimental ones. This is interpreted as evidence 
obtained from antiprotonic atoms for too rapid a decrease of these theoretical nucleon densities as a function of the 
radial distance due to a too small diffuseness of these densities. A similar conclusion was previously obtained from 
the analysis of nucleon elastic scattering. 
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Tables 



TABLE I: Target properties and number of antiprotons used. 



Target 


thickness d 


enrichment 


number of p (10 ) 




(mg/cm 2 ) 


(%) 




2U8 Pb 


130.4 


99.1 


17 


209 Bi 


132.7 


nat. 


1.4 
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TABLE II: Measured relative X-ray intensities, normalized to the tran- 
sition n= 13 — > 12 (average value from three detectors). 



Transition 




2,19 Bi 


10-> 


9 


69.0 ± 2.3 


63.6 ± 3.3 




10 


99.6 ± 3.2 


97.8 ± 3.6 


12-> 


11 


103.8 ± 4.1 


102.1 ± 4.7 


lo — ► 


1 9 


100.0 ± 4.6 


100.0 ± 4.9 


14-> 


13 


93 5 ± 5 


93 8 ± 5 4 


15-> 


14 


81.3 ± 3.7 


81.2 ± 4.0 


16-> 


15 


59.4 ± 2.7 


61.4 ± 3.1 


17-> 


16 


54.8 ± 10.9* 


80.3 ± 10.4* 


18-> 


17 


73.4 ± 3.6 


71.2 ± 4.0 


19-> 


18 


56.0 ± 2.6 


56.6 ± 2.9 


20^ 


19 


48.2 ± 2.3 


53.4 ± 4.8 


21-> 


20 


44.9 ± 3.3 


67.8 ± 6.3 


ll-> 


9 


3.3 ± 0.3 




12-f 


10 


7.6 ± 0.5 


5.3 ± 0.8 


13-f 


11 


10.1 ± 0.3 


9.1 ± 0.5 


14-> 


12 


11.2 ± 0.6 


10.4 ± 0.8 


15-> 


13 


10.6 ± 0.5 


11.2 ± 0.8 


16-> 


14 


10.5 ± 0.5 


10.7 ± 0.8 




1 ^ 

LO 


12.2 ± 0.6 


11.4 ± 0.8 


18-> 


16 


12.0 ± 0.4 


11.2 ± 0.6 


19-> 


17 


11.8 ± 0.6 


11.0 ± 0.8 


20^ 


18 


10.1 ± 0.5 


11.0 ± 0.8 


21-> 


19 


5.2 ± 2.8* 


1.3 ± 1.0 * 


22^ 


20 


0.0 ± 7.9* 


0.0 ± 7.8* 


23-> 


21 


4.2 ± 0.2 


9.6 ± 1.1 


24^ 


22 


3.5 ± 0.3 


4.4 ± 0.6 


25^ 


23 


5.0 ± 0.3 


5.1 ± 0.7 


14-> 


11 


3.1 ± 0.3 




15-> 


12 


2.6 ± 0.2 


3.3 ± 1.1 


16-f 


13 


2.4 ± 0.2 


3.0 ± 0.6 


17-f 


14 


4.2 ± 0.3 


2.4 ± 0.5 


18-c 


15 


3.7 ± 0.2 


3.1 ± 0.5 


19-f 


16 


3.5 ± 0.2 


2.8 ± 0.5 


on ». 


1 7 


3.5 ± 0.2 


2.8 ± 0.4 


21-f 


18 


3.7 ± 0.2 


2.9 ± 0.3 


22^ 


19 


3.9 ± 0.2 


3.0 ± 0.4 


23-> 


20 


3.3 ± 0.2 


2.4 ± 0.5 


24^ 


21 


1.4 ± 0.3 


3.4 ± 0.7 


25^ 


22 


3.0 ± 0.2 


2.0 ± 0.4 


26^ 


23 


10.3 ± 0.6 


11.5 ± 0.8 


27^ 


24 


9.8 ± 0.5 


3.5 ± 0.7 


28^ 


25 


3.5 ± 0.3 




17-» 


13 


3.3 ± 0.3 




18-> 


14 


1.2 ± 0.2 




19-> 


15 


1.6 ± 0.2 




20^ 


16 


1.7 ± 0.2 




21-> 


17 


2.1 ± 0.2 




22^ 


18 


2.8 ± 0.2 




23^ 


19 


1.2 ± 0.1 




24^ 


20 


1.3 ± 0.2 




25-> 


21 


1.1 ± 0.2 





* admixtures of electronic X-rays from the same atom 
and from the (Z±l) atoms were subtracted. 
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TABLE III: Measured shifts of the (n, I — 10, 9) (e u ) and (n, I = 9, 8) (e;) levels in the antiprotonk 


208 Pb and 209 Bi atoms. 


Target e+ (eV) e" (eV) e+ (eV) 


ef ( eV ) 


2UB Pb 34 ± 16 28 ± 17 102 ± 28 
209 Bi 23 ±20 -8 ± 21 29 ± 72 


73 ± 29 
-2 ± 73 


TABLE IV: Measured absorption widths of the fine structure components of the (n, Z = 9,8) level 
and 209 Bi atoms. 


in the antiprotonic 208 Pb 


Target T+ (eV) 


rr (eV) 


2U8 Pb 320 ± 35 
209 Bi 557 ± 68 


302 ± 38 
448 ± 74 


TABLE V: Radiation width T em and Auger width T Auger for the n=10 levels, where the strong 
determined via the intensity balance. 


interaction width T u was 


Target Y em (eV) V Auger (eV) T+ (eV) 


I 1 " (eV) 


/U8 Pb 12.59 0.139 5.3 ± 1.1 
209 Bi 13.27 0.141 6.1 ± 1.7 


6.6 ± 1.3 
7.8 ± 1.9 



TABLE VI: Comparison of 2pF proton and neutron distributions in 208 Pb (all parameters in fm). 





protons 


neutrons 








a p 


Op 




a n 


c„ 






ATnp 




(a) 








(a) 


(b) 










(a) 


(ft) 


SkP 


5.465 


5.489 


0.437 


6.768 


5.610 


5.625 


0.537 


6.789 


0.100 


0.021 


0.145 


0.136 


SkX 


5.441 


5.443 


0.424 


6.726 


5.597 


5.597 


0.510 


6.799 


0.086 


0.073 


0.156 


0.154 


DD-ME2 


5.460 


5.472 


0.444 


6.736 


5.653 


5.657 


0.561 


6.789 


0.117 


0.053 


0.193 


0.185 


Fricke (c) 




5.436 


0.446 


6.684 


















Experiment 1 ^ 












5.596 


0.571 


6.684 


0.125 


0.0 (e) 




0.16(2) 



( a > calculated from the theoretical distributions; 
' ' calculated from fit parameters: a p ,c p , a n , c n ; 

( c ) point proton values obtain from the Fricke [45] charge distribution using Oset's [481 transformation 
formulae; 

^ from 2pF fit to the experimental width using the Batty potential [121 : 
assumed. 



TABLE VII: Parameters of the annihilation probability distribution calculated using DD-ME2 density and the zero-range Batty 
potential. 



Distribution 


Level 


FWHM 


radial parameter (fm) 
most probable median 


average 


DD-ME2 


up 


1.5 


8.3 8.7 


9.1 




low 


1.2 


8.5 8.6 


8.8 
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TABLE VIII: Comparison of experimental level widths and shift in Pb with those calculated using theoretical neutron and 
proton distributions. 





~^low 


Tup 


e 




eV 


eV 


eV 


Experiment 


312(26) 


5.9(8) 


88(20) 


Batty potential 








SkP 


274 


5.2 


14 


SkX 


231 


4.2 


16 


DD-ME2 


315 


6.2 


12 


Friedman potential 








SkP 


278 


5.3 


6 


SkX 


244 


4.5 


7 


DD-ME2 


307 


6.1 


2 
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FIG. 1: Summary of shifts and widths measured for Pb. 



February 8, 2008 



14 



Energy (keV) 

50 100 150 200 250 300 350 400 450 



50000 



g 

U 40000 
o 
<u 

Oh 



o 
U 



30000 



20000 - 



10000 



18->17 

19— >18 
20^19 

V 1 



i i i r 



17^16 



15— >14 



J L 



14— >13 
(18-»16) 
13->12 



12— >11 



11-»10 




10^9 
1,(13^11) 



250 500 750 1000 1250 1500 1750 2000 2250 2500 
Channel number 

FIG. 2: Antiprotonic X-ray spectrum from 208 Pb measured with the HPGe detector of 19% relative efficiency. 
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FIG. 3: Part of the antiprotonic X-ray spectrum measured for 208 Pb using the detector with the 1035 mm 2 x 14 mm crystal. 
The fit to the broadened 10— »9 transition is also shown. The 13— >11 line is admixed to the 10^9 line. 
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FIG. 5: Width and shift of the levels (n, / = 9, 8) and widths of the levels (n, I = 10, 9) plotted versus Z, for 208 Pb and 209 Bi. 
Positive level shifts correspond to repulsive interactions. The presented lower level shift and width for 209 Bi are not corrected 
for the hyperfine contribution. 



February 8, 2008 



16 



in 

as 

M 

E 



Q. 
Q. 




FIG. 6: Comparison of 208 Pb charge-density distributions deduced from various compilations. The plotted charge density 
distributions are normalized to the o ne g iven by Fricke I45|. Other charge distributions (charge rms radii in parenthesis) are 
taken from: (A) gl_[5.521 fm), (B) HI (5.515 fm), (C)[63 (5.510 fm), (D) (5.503 fm), (E) [«| (5.46 fm). The charge 
rms radius given in [45| is 5.504 fm. 
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FIG. 7: Left panel: two-parameter Fermi (2pF) distribution fitted in the range 1-10 fm to proton and neutron distributions 
calculated using the relativistic mean-field theory (RMF) with DD-ME2 parametrization 154]- Points - calculated distributions, 
continuous line - 2pF fit. Right panel: relative differences of fitted (2pF) and calculated (DD-ME2) densities. Continuous line 
- protons, dashed line - neutrons. 
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FIG. 8: Sum of the theoretical neutron and proton densities for the radial distances at which the antiproton annihilation in 
208 Pb nucleus is significant. (A) - HF model with SkX parameters, (B) - HFB model with SkP parameters, (C) - RMF model 
with DD-ME2 parameters. 
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FIG. 9: Normalized neutron to proton point-nucleon density ratios as a function of the distance from the nuclear center in 
208 Pb. Left panel (cf. Section HV C|) : (A) - two-parameter Fermi distribution with Ar np =0.16 fm and Ac np =0, (B) - SkP 
parametrization, (C) - SkX parametrization, (D) - DD-ME2 parametrization. Right panel (cf. Section [IV G[) : density ratios 
deduced from two-parameter Fermi distributions: (A) - the same as on the left panel, (E) - Ar np =0.16 fm, Ac np =0.1 fm, (F) 
- Ar np =0.16 fm, Ac np =0.2 fm, (G) Ar n p=0.17 fm, Ac np =0.13 fm (cf. Sec. IIV G 3[l . The cross at 8 fm represents the halo 
factor from Bugg's experiment [43], the cross at 9.1 fm is the interpolated halo factor (see Sec. IIV F)l . 
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FIG. 10: The p annihilation probability (arbitrary units) from the upper state (dashed line) and the lower state (continuous line) 
for the antiprotonic 208 Pb atom. The matter density (DD-ME2) in 208 Pb is also shown. The calculations of the annihilation 
probability were done with the Batty zero-range potential and 2pF parametrization of the DD-ME2 density (cf. Table IVT)) . 
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FIG. 11: A straight line fitted to the experimentally determined halo factors [1,0, HH], plotted as a function of the target-nucleus 
neutron binding energy. For 208 Pb the neutron binding energy is equal to 7.4 MeV and the interpolated halo factor is 2.8 ±0.4. 
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FIG. 12: The relation between the difference of the diffuseness parameters a and half-density radii c of neutron and proton 
distributions in 208 Pb deduced from the experimental antiprotonic level widths using the zero-range interaction Batty potential 
(Ref. [l^])- The experimental uncertainties of the level widths are reflected by an almost constant Aa np uncertainty (not shown) 
of about ±0.015 fm. 
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FIG. 13: The relation between the difference of the neutron and proton half-density radii Ac„ p and the deduced difference 
between rms radii of neutron and proton distributions Ar np of 208 Pb for the zero-range interaction Batty potential. The dashed 
line and shaded region shows the weighted average and the error, respectively, of Ar np determined in the hadron scattering 
experiments (see Ref. [H). 
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FIG. 14: The same as Fig. 1121 but for the finite-range Friedman potential 
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